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Abstract. It is known that the famous Feigenbaum-Coullet-Tresser period 
doubling universality has a counterpart for area-preserving maps of IR 2 . A 
rcnormalization approach has been used in (Eckmann et al 1982) and (Eckmann et 
al 1984) in a computer-assisted proof of existence of a "universal" area-preserving 
map F* — a map with orbits of all binary periods 2 k , k £ N. In this paper, we 
consider maps in some neighbourhood of F* and study their dynamics. 

We first demonstrate that the map F* admits a "bi-infinite heteroclinic 
tangle": a sequence of periodic points {z k }, k £ Z, 

\z k \ fc -^5° 0, \z k \ oo, (1) 

whose stable and unstable manifolds intersect transversally; and, for any N £ N, 
a compact invariant set on which F* is homeomorphic to a topological Markov 
chain on the space of all two-sided sequences composed of N symbols. A corollary 
of these results is the existence of unbounded and oscillating orbits. 

We also show that the third iterate for all maps close to _F* admits a horseshoe. 
We use distortion tools to provide rigorous bounds on the Hausdorff dimension of 
the associated locally maximal invariant hyperbolic set: 

0.7673 > dim H (C F ) >e^ 0.00013 e^ 7499 . 
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1. Introduction 

Following the pioneering discovery of the Feigcnbaum-Coullet-Tresser period doubling 
universality in unimodal maps (Feigcnbaum 1978), (Feigenbaum 1979), (Tresser and 
Coullet 1978) universality has been demonstrated to be a rather generic phenomenon 
in dynamics. 

To prove universality one usually introduces a renormalization operator on a 
functional space, and demonstrates that this operator has a hyperbolic fixed point. 

Such renormalization approach to universality has been very successful in one- 
dimensional dynamics, and has led to explanation of universality in unimodal maps 
(Epstein 1986), (Epstein 1989),(Lyubich 1999), critical circle maps (de Faria 1992, de 
Faria 1999), (Yampolsky 2002), (Yampolsky 2003) and holomorphic maps with a Siegel 
disk (McMullen 1998), (Yampolsky 2007), (?). 

Universality has been abundantly observed in higher dimensions, in particular, in 
two and more dimensional dissipative systems (cf. (Collet et al 1980), (Sparrow 1982)), 
in area-preserving maps, both as the period-doubling universality (Derrida and 
Pomcau 1980), (Hcllcman 1980), (Bcncttin et al 1980), (Collet et al 1981), (Eckmann 
et al 1982), (Eckmann et al 1984), (Gaidashev and Koch 2008) and as the universality 
associated with the break-up of invariant surfaces (Shcnkcr and Kadanoff 1982), 
(MacKay 1983), (MacKay 1982), (Mchr and Escande 1984), and in Hamiltonian flows 
(Escande and Doveil 1981), (Abad et al 2000), (Abad et al 1998), (Koch 2002), (Koch 
2004), (Koch 2004), (Gaidashev and Koch 2004), (Gaidashev 2005), (Kocic 2005) . 

It has been established that the universal behaviour in dissipative and 
conservative higher dimensional systems is fundamentally different. The case of the 
dissipative systems is often reducible to the one-dimensional Feigcnbaum-Coullet- 
Tresser universality ((Collet et al 1980), (de Carvalho et al 2005), (Lyubich and 
Martens 2008)). Universality for highly dissipative Henon like maps 

F(x,u) = (f(x) - e(x,u),x), 

where / is a unimodal map, e - sufficiently small together with its derivatives, has 
been demonstrated in (de Carvalho et al 2005). Specifically, it has been shown that 
these maps arc in the Feigcnbaum-Coullet-Tresser universality class, and that the 
degenerate map 

F*(x,u) = (f*(x),x), 

where /* is the Feigcnbaum-Coullet-Tresser fixed point, is a renormalization fixed 
point in that class. The authors of (de Carvalho et al 2005) have also constructed a 
Cantor set (an attractor) for infinitely rcnormalizablc maps, and demonstrated that 
for such Hcnon-likc maps universality coexists with non-rigidity: the Holder exponent 
of the conjugacy between the actions of any two infinitely renormalizablc maps with 
nonequal average Jacobians on their Cantor sets has an upper bound less than 1. 

The case of area-preserving maps seems to be very different, and at present 
there is no deep understanding of universality in conservative systems, other than 
in the "trivial" case of the universality for systems "near integrability" (Koch 2002), 
(Koch 2004), (Gaidashev 2005), (Kocic 2005), (Khanin et al 2007). 

An infinite period-doubling cascade in families of area-preserving maps was 
observed by several authors in early 80's (Derrida and Pomeau 1980), (Helleman 1980), 
(Benettin et al 1980), (Bountis 1981), (Collet et al 1981). The period-doubling 
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phenomenon can be illustrated with the area-preserving Hcnon family (cf. (Bountis 
1981)) : 

H a (x, u) — (—it + 1 — ax 2 , x). 



Maps H a have a fixed point ((—1 + yl + a)/a, (—1 + \/TTa)/a) which is stable 
for — 1 < a < 3. When ai = 3 this fixed point becomes unstable, at the same time 
an orbit of period two is born with H a (x±,x T ) — (x T ,x±), x± = (1 ± \J a — 3)/a. 
This orbit, in turn, becomes unstable at ai = 4, giving birth to a period 4 stable 
orbit. Generally, there exists a sequence of parameter values a&, at which the orbit of 
period 2 fc_1 turns unstable, while at the same time a stable orbit of period 2 k is born. 
The parameter values accumulate on some . The crucial observation is that the 
accumulation rate 



lim 

k — *oo ak-\-l 



Ofc ^ Ctk-1 



= 8.721. 



Ofc 



(2) 



is universal for a large class of families, not necessarily Henon. 

Furthermore, the 2 h periodic orbits scale asymptotically with two scaling 
parameters 

A =-0.249..., M = 0.061... (3) 

To explain how orbits scale with A and li we will follow (Bountis 1981). Consider 
an interval (ak,ak+i) of parameter values in a "typical" family F a . For any value 
a 6 (afc, dk+i) the map F a possesses a stable periodic orbit of period 2 k . We fix some 
a/c within the interval (ak,a,k+i) in some consistent way; for instance, by requiring 
that the restriction of F 2 k to a neighbourhood of a stable periodic point in the 2 fc - 
periodic orbit is conjugate, via a diffeomorphism H^, to a rotation with some fixed 
rotation number r. Let p' k be some unstable periodic point in the 2 fc ~ 1 -periodic orbit, 
and let pk be the further of the two stable 2 fe -periodic points that bifurcated from p' k . 
Denote with dk = \p' k — Pk | > the distance between p^ and p' k . The new elliptic point 
Pk is surrounded by invariant ellipses; let Ck be the distance between pk and p' k in the 
direction of the minor semi-axis of an invariant ellipse surrounding p k , see Figure [T] 
Then, 



— = — lim 
A k 



dk 
00 d k+ i 



A 



Pk 



= — lim 

[X k^oo p k+1 



1 ,. c k 

Y2 = , llm • 

A fc— >oo Ck+1 



where pk is the ratio of the smaller and larger eigenvalues of DHk{pk)- 




c k 



Figure 1. The geometry of the period doubling, pj. is the further elliptic point 
that has bifurcated from the hyperbolic point p' k . 
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This universality can be explained rigorously if one shows that the renormalization 
operator 

R[F]=Kp 1 oFoFoK F , (4) 

where is some independent coordinate transformation, has a fixed point, and the 
derivative of this operator is hyperbolic at this fixed point. 

It has been argued in (Collet et al 1981) that Ap is a diagonal linear 
transformation. Furthermore, such Ap has been used in (Eckmann et al 1982) 
and (Eckmann et al 1984) in a computer assisted proof of existence of a reversible 
renormalization fixed point and hyperbolicity of the operator R. 

An exploration of a possible analytic machinery has been undertaken in 
(Gaidashcv and Koch 2008) where it has been demonstrated that the fixed point 
is very close, in some appropriate sense, to an area-preserving Henon-like map 

H*(x, u) — (4>(x) — u, x — 4>{4>{x) — u)), (5) 

where cf) solves the following one-dimensional problem of non-Feigenbaum type: 

<j>{y) = \<i>{4>^y))-y- (6) 

In this paper we will study the dynamics of the renormalization fixed point F* and 
maps in some neighbourhood of F„ . We would like to emphasize that the preservation 
of area for such maps leads to one important difference between the case at hand 
and the highly dissipative case of Henon-like maps: the former maps have no obvious 
invariant subsets in their domain. In particular, this means, that the construction of 
invariant Cantor sets carried out in (de Carvalho et al 2005), which is using existence 
of invariant subsets in a crucial way, is not readily applicable to the case of area- 
preserving maps. 

To construct hyperbolic sets, we will use the idea of covering relations (see, 
e.g. (Zgliczyhski and Gidea 2004, Zgliczyhski 2009)) in rigorous computations. The 
Hausdorff dimension of the hyperbolic sets will be estimated with the help of the 
Duarte Distortion Theorem (see, e.g. (Duarte 2000)) which enables one to use the 
distortion of a Cantor set to ultimately find bounds on the dimension. 

The structure of the paper is as follows: we begin by recalling the basic properties 
of area-preserving reversible maps in Section 2. In Section 3 we introduce some 
notation and recall some standard definitions from hyperbolic dynamics. In Section 4 
we study the domain of analyticity of maps in a neighbourhood of the renormalization 
fixed point, and analytic continuation of the renormalization fixed point. Section 5 
consists of the statements of our main theorems. In Section 6, we recall the definition 
and the main properties of the covering relations. In Section 7 we prove that any area- 
preserving reversible map in a neighbourhood of the renormalization fixed point has a 
transversal homoclinic orbit in its domain of analyticity. In Section 8 we construct a 
heteroclinic tangle for the renormalization fixed point. From its existence the existence 
of unbounded, and oscillating trajectories follow. In Section 9 we recall the Duarte 
Distortion Theorem. In Sections 10 and 11 we pr ove that the third iterate of any 
area-preserving reversible map in a neighbourhood of the renormalization fixed point 
has a horseshoe, and compute bounds on its Hausdorff dimension using the Duarte 
distortion theorem. 

In a satellite paper (Gaidashev and Johnson 2009), we prove that infinitely 
renormalizable maps in the neighbourhood of existence of the hyperbolic set for the 
third iterate also admit a "stable" set. This set is a bounded invariant set, such that 
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the maximal Lyapunov exponent for the third iterate is zero. In (Gaidashev and 
Johnson 2009) we provide an upper bound on the Hausdorff dimension of the stable 
set, and prove that the Hausdorff dimension is constant for all maps in some subset 
of infinitely renormalizable maps. 



2. Renormalization for area-preserving reversible maps 

An "area-preserving map" will mean an exact symplectic diffeomorphism of a subset 
of K 2 onto its image. 

Recall, that an area-preserving map can be uniquely specified by its generating 
function S: 

\ F / 
x i i > I V 



-Si(x,y)J \S 2 (x,y)J' St dlS ' ^ 
Furthermore, we will assume that F is reversible, that is 

T o F oT = F~ x , where T(x, u) = (x, -u). (8) 
For such maps it follows from Q that 

Si(y, x) = S 2 (x, y) = s(x, y), (9) 

and 

- s (y,x)J \s{x,y)J 
It is this "little" s that will be referred to below as "the generating function". 
It follows from (JoJ) that si is symmetric. If the equation —s(y,x) — u has a unique 
diffcrcntiablc solution y = y{x, u), then the derivative of such a map F is given by the 
following formula: 



DF(x,u) 



s 2 (y(x,u),x) 



s 1 (y(x,u),x) s 1 (y(x,u),x) 

Sl (x, y(x, u)) - s 2 (x, y(x, u)) s ^ u \ x \ _ Mxv(xu)) 

1 ^ ?y V ' '/ z ^ ' y ^ ' // 8i(y(x,u),x) si{y(x,u),x) 



(ii) 



We will now derive an equation for the generating function of the renormalizcd 
map Ap o F o F o . 

Applying a reversible F twice we get 

x ' \~ ( z ' \ = ( z ' \~ ( y ' 

-s(z',x')J \s(x',z')J {-s(y',z')J \s(z',y>) 
It has been argued in (Collet et al 1981) that 

Ap(x,u) — (\x,[iu). 
We therefore set (x' , y') = (Xx, Xy), z'(Xx, Xy) = z(x, y) to obtain: 

* \ Af ( \ FoF f Xy f y \ 

l-i ■ \ l ^ ^ (12) 



■±s(z,\x)J \-s(z,Xx)J \s(z,Xy)J \is(z,Ay) 



where z{x, y) solves 

s(Xx, z(x, y)) + s(Xy, z(x, y)) = 0. (13) 
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If the solution of (13) is unique, then z(x,y) — z(y,x), and it follows from (12 1 
that the generating function of the renormalized F is given by 

s(x,y) = ^ 1 s(z(x,y), Xy). (14) 

One can fix a set of normalization conditions for s and z which serve to determine 
scalings A and p as functions of s. For example, the normalization 

s(l,0) = 

is reproduced for s as long as 

2(1,0) =2(0,1) = 1. 
In particular, this implies that 

s(A,l) + s(0,l) = 0. 

Furthermore, the condition 

0is(l,Q) = l (15) 
is reproduced as long as 

/i = 012(1,0). 

We will now summarize the above discussion in the following definition of the 
renormalization operator acting on generating functions originally due to the authors 
of (Eckmann et al 1982) and (Eckmann et al 1984): 



Definition 2.1 



K EKW [s)(x,y) = v 1 s(z(x,y),Xy), (16) 

where 

= s(Xx, z(x, y)) + s(Xy, z(x, y)), (17) 

= s(A,l) + s(0, 1) and /i = <9iz(l,0). (18) 



Definition 2.2 The Banach space of functions s(x,y) = 53ij=o c ij xl U^ > analytic on 
a bi-disk 

|x-0.5| <p, |y-0.5| <p, 

for which the norm 

00 

\\s\\ P = E \^\p l+J 

ij=0 

is finite, will be referred to as A(p). A. s (p) will denote its symmetric subspace 
{s e A(p) : si(x,y) = si(y,x)}. 



As we have already mentioned, the following has been proved with the help of a 
computer in (Eckmann et al 1982) and (Eckmann et al 1984): 



Dynamics of the Universal Area- Preserving Map: Hyperbolic Sets 6 

Theorem 2.3 There exist a polynomial s a G A s (p) and a ball B r (s & ) C A s {p), 
r = 6.0 x 1CP 7 . p — 1.6, such that the operator IZekw * s well-defined, analytic and 
compact on B r . 

Furthermore, its derivative T)TLekw\b t has exactly two eigenvalues 8± and 62 of 
modulus larger than 1, while 

spec(DTZ E Kw\B r ) \ {Si,S 2 } C {z e C : \z\ < v < 1}. 

Finally, there is an s* G B r such that 

Ti-EKw[s*] = s*. 

The scalings A* and /i* corresponding to the fixed point s* satisfy 

A, e [-0.24887681,-0.24887376], (19) 
£ [0.061107811,0.061112465]. (20) 



Remark 2.4 The radius of the contracting part of the spectrum spec(DlZEK\v( s *) \ 
{Si, 62} has been estimated in (Eckmann et al 1984) to be is = 0.8. 

It follows from the above theorem that there exist a codimension 2 local stable 
manifold Wf oc (s*) C B r . 



Definition 2.5 A reversible map F of the form (10) such that s £ W* oc (s*) is called 
infinitely renormalizable. The set of all reversible infinitely renormalizable maps is 
denoted by W. 



3. Some notation and definitions 



We will use the following notation for the sup norm of a function h and a 
transformation H defined on some set 5cl 2 or C 2 : 

\h\ s = sup {\h\}, (21) 
|tf|s = max{ sup {\V X H\}, sup {\V U H\}}, (22) 

where V x and V u are projections on the corresponding components. 
We will also use the notation | • | for the I2 norm for vectors in K 2 . 
The interval enclosures of A* and /i* will be denoted 

A* = [A_, A + ]; A_ = -0.24887681, A+ = -0.24887376, (23) 

/a, = M- = 0.061107811, ^+ = 0.061112465. (24) 

The corresponding interval enclosure for the linear map A* will be denoted A*; 
if (x, u) e C 2 , then 

A* (x, u) = { ( \x, fiu) eC 2 :AeA„(i£fi,}. (25) 

The bound on the fixed point generating function s* will be called s*: 

s* = {s e A s {p) : \\s - s & \\ p < r = 6.0 x lO" 7 } , (26) 

while the bound on the renormalization fixed point F* will be referred to as F*: 

F, = {F : (x, -s(y, x)) ~ (y, s(x, y)) :s£s'}, (27) 
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where s a is as in Theorem |2.3| the third iterate of this bound will be referred to as 
G*. With DG* : p i— > [DG»(p)] we denote an interval matrix valued function such 
that 

[DG(p)]ij e [DG.ipflij, for all GgG„ P £ V 3 , 

where T>% is the domain of G*, and the bound on the operator norm of DG for 
G = F o F o F, F e on a set S will be denoted 

llDG.Hs = sup {\\D{F o F o F)\\ s } . 
fgf, 

Given a non-empty open set T> C C™ we will denote by Or(D) the set of reversible 
area-preserving maps F :T> i— > C™, analytic on T>. 

We will proceed with a collection of classical notations (see, eg, (Katok and 
Hasselblatt 1995)) relevant to our following discussion. 

Definition 3.1 (Hyperbolic set) Let M. be a smooth manifold, and let F be a 
diffeomorphism of an open subset IA C M. onto its image. 

A set C is called hyperbolic for the map F if there is a Riemannian metric on a 
neighbourhood 11 of C , and < 1 < 8, such that for any p € C and n G N the tangent 
space Tpn^lA admits a decomposition in two invariant subspaces: 

T Fn[p) M = E+® E-, DF(F n (p))E± = E± +v 

on which the sequence of differentials is hyperbolic: 

\\DF{F n {p))\ E -\\<(3, \\DF-\F n {p))\ E+n J<5-\ 

Definition 3.2 (Locally maximal hyperbolic set) LetC be a hyperbolic set for F : hi 
M.. Lf there is a neighbourhood V of C such that C = n ne zF n (V), then C is called 
locally maximal. 

Definition 3.3 (Bernoulli shift) Let {0,1, . . . , N — 1} Z be the space of all two-sided 
sequences of N symbols: 

{0, 1, . . . ,N - if = {oj = (. . .,w_i,Wo,wi, ...)■• Wi € {0,1,... ,N- l},i e Z}, 

Define the Bernoulli shift on {0, 1, . . . , N — 1} Z as 

v n (uj)=J, oj' n = uj n+ i. 

Definition 3.4 (Topological Markov chain) Let A = (a>i,j)ij^o be an N x N matrix 
whose entries are either or 1. Let 

{0,1, . . . ,N - 1}\ = {w £ {0, 1, . . . , N - 1} Z : a Un . Un+1 = l,n € Z}. 

The restriction 

(J iv|{0,l,...,iV-l}Z = a A 

is called a topological Markov chain determined by A. 
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Definition 3.5 (Homoclinic and heteroclinic points) Let F : X i— ► X be a 

homeomorphism on a metric space (X,d). A point x £ X is said to be homoclinic 
to the point y G X if 

lim d(F n (x),F n (y)) = 0. 

| n| — >oo 

A point x is said to be heteroclinic to points yi and y 2 if 

lim d(F n {x),F n (y 1 )) = lim d(F n (x), F n (y 2 )) = 0. 



n — > — oo 



In this case we will also say that there exists a heteroclinic orbit between points y\ and 

If A4 is a differentiable manifold, F £ Diff (M) and x £ M. is a hyperbolic fixed 
point of F , then we say that q £ M is a transversal homoclinic point to x if it is a 
point of transversal intersection of the stable and unstable manifolds of x. 

Definition 3.6 Let X be a metric space. If S C X , and d £ [0, oo), the d-dimensional 
Hausdorff content of S is defined as 



Cf(S) = inf |^]diam(5 l ) d : {S, } is a cover of S j . 



(28) 



The Hausdorff dimension of S is defined as 



dim ff (S) = inf {d > : Cf{S) = 0} . (29) 



4. Domain of analyticity of and G* 

F is defined implicitly by the generating function s, its domain is given as in (Eckmann 
et al 1984): 

V s = {(x, y) e C 2 : | a; - 0.5| < 1.6, \y - 0.5| < 1.6}. (30) 

To find the domain of F € F*, we note that its second argument is equal to —s(y, x), 
for some s £ s* (see (10)). Thus, the domain of F, T), is given by: 

V = {(x,u) £ C 2 : u = -s(y,x), s £ s*(y,x), \x - 0.5| < 1.6, \y-0.5\ < 1.6}. (31) 

We denote by 

T> = {(x,u) £V :3a; = 3u = 0}, (32) 
the real slice of T>. 

To solve nonlinear equations on the computer we use the interval Newton 
operator, see e.g. (Neumaier 1990). 

Definition 4.1 (Interval Newton Operator ) Let F : V C E™ ->■ R™, and let DF be 

an interval matrix valued function such that [DF(x)]ij £ [D(F)(x)]ij, for all x £T>. 
Let X C T> be a Cartesian product of finite intervals, x £ X, and assume that if 
A £ DF(X), then A is non-singular. We define the interval Newton operator as: 

N(F,X,x) = x-{DF)- 1 (X)F(x). 

The main properties of N is that if N(F, X, x) C intX, then there exists a unique 
solution to F{x) =0mX, which is contained in N(F,~X., x) , and if N(F, X, i)nX = 0. 
then there is no solution to F(x) =0 in X. 
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Lemma 4.2 There exists a non-empty open set T> , 

V X V c {i £ I : |ar — O.Sj < 1.6}, 

such that for every (x,u) E T> there exists a unique solution of the interval Newton 
operator of the function h( x>u -\(y) = u + s(y,x), s € s*(y,x), that satisfies \y — 0.5| < 
1.6. 

Proof. Set Y = {y e K : \y - 0.5j < 1.6}. Given (x,u), let N(h {XtU) , Y, y) be the 
interval Newton operator (for some appropriately chosen y S Y). We have verified 
that there exists a non-empty set T>, such that for all (x, it) G I? 

N(h {x>u) ,Y,y) CY. 

This verification is implemented in the program f indDomain of (Progsl 2009) It fol- 
lows, see e.g. (Neumaier 1990), that there is a unique y, such that (x,y) is in the real 



slice of V s and h( x ^ u \(y) = 0. Thus, y is defined as a function of (x, u), and by (31 1 
(x,u) € V. 

□ 

Clearly, D C V. 

The generating function s is analytic on a bi-disk, preferably F should have a 
similar property. From the above construction, we can at least show that F is defined 



on a complex neighbourhood of T> from the above Lemma 4.2 



Lemma 4.3 T> contains an open complex neighbourhood of the set T>. 

Proof. N(h( x u } , y) is only well defined if s* i(y,x) ^ 0, so that the condition of the im- 
plicit function theorem is automatically satisfied on the solution set of N(hr x >u \, Y, y). 
Since s*i is analytic on the bi-disk {|x — 0.5| < 1.6} x {\y — 0.5 1 < 1.6}, there is an 
open neighbourhood in C 2 of the solution set where s*i(y, x) ^ 0. It follows from the 
implicit function theorem that is analytic on this neighbourhood. 

□ 

Generally, we will denote the third iterate of a map F £ Or(T>) as G. The domain 
of G, 2?3, is given by 

V 3 = F~ 2 (V"nV), (33) 

where 

V" = F(V' n V) , T>' = F(D). 

Using the program f indDomain, we have verified that the real slice Th of 2?3 is 
an open non-empty set. Approximations of T> and T) 3 are shown in Figured] We note 
that by using the renormalization equation F* = R[F*], F* has, for any k, an analytic 
continuation to domains 

v k = A- 1 (F- 1 (F„(v k - 1 ) n X> fe_1 )), V° = V, 

while G* has an analytic continuation to 

V^A-^G-^Vl-^nVt 1 )), -D°s = ^ 3 . 

The real slices of the domains T>® and T>\ are given in the Figure [3] 
We will conclude this Section with the following 
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2.5 



Figure 2. The real slices of the domains of F„ (blue) G» (red). 



Lemma 4.4 All maps FgF, posses a hyperbolic fixed point pq = po(F), such that 

1) VxPa e (0.57761843,0.57761989), and V u p = 0; 

2) DF(po) has two real eigenvalues: 

e+ g (-2.05763559,-2.05759928), 
e_ G (-0.48601715,-0.48598084). 

Proof. The bound on the fixed point has been obtained with the help of the interval 
Newton method. Hyperbolicity has been demonstrated by computing a bound on DF 



using formula (11) 



□ 



5. Statement of the main results 

We will now summarize our main findings. 
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Figure 3. The real slices of the domains Dg (red) and (blue). 



Our first theorem describes the two-sided heteroclinic tangle for the fixed point 
map F*. 

Main Theorem 1 The renormalization fixed point F* has the following properties: 

1) F* possesses a point which is transversally homoclinic to the fixed point po; 

2) there exists a positive integer n such that for any negative integer k the map i 7 "™ 
has a heteroclinic orbit Ok between the periodic points A^(po) and A^ +1 (po), and 
for any positive integer k the map F™' 2 has a heteroclinic orbit Ok between the 
periodic points Aj(po) and A^~ 1 (pq); 

3) for any N G N and MfN there exists an integer n and an invariant set Zp™, 
such that ^ 

F"\z pn , ^a\{-n,-n+i,—,m-i,m} z > 
* homeo 

where cfa is the topological Markov chain defined by a (N + M + l) x (N+M + l) 
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tridiagonal matrix 

' 





1110 
111 
11 

4) for any R > e > there exists a point p eT> and an n G N such that \F™(p) | > R 
fan unbounded orbit ), and a point q G V and n,m G N such that \F™(q)\ > R 
and \F™(q) \ < e (an oscillating orbit). 

Our second result demonstrates that all locally infinitely renormalizable maps 
admit a hyperbolic set in their domain of analyticity. 

Main Theorem 2 Any FeF, admits a hyperbolic set C G C V 3 for G = F o F o F; 

whose Hausdorff dimension satisfies: 

0.7673 > dim H (C G ) > e, 
where e ~ 0.00013 e -7499 is strictly positive. 

6. Topological tools 

The main tools of our proofs are covering relations (Zgliczyhski 1997, Zgliczyhski and 
Gidea 2004) and cone conditions (Kokubu et al 2007, Zgliczyhski 2009), see also (Galias 
and Zgliczyhski 2001) for proving the existence of homoclinic and hctcroclinic orbits. 
To make the present paper reasonably self-contained we include a brief introduction 
to the necessary concepts. 

6.1. H-sets and Covering relations 

The notion of an h-set and a covering relation hrst appeared in (Zgliczyhski 1997), 
the most thorough treatment is (Zgliczyhski and Gidea 2004). The basic idea is to 
construct computable conditions for the existence of a scmi-conjugacy to symbolic 
dynamics. This is done by constructing h-sets, i.e. hyperbolic-like sets, that cross 
each other in a (topologically) non-trivial way. We denote by B n (c,r), the open ball 
in R™ with centre c and radius r, and S n (c,r) = dB n+ i(c,r). 

Definition 6.1 An h-set is a quadruple consisting of 

• a compact subset \N\ o/R™, 

• a pair of numbers u(N),s(N) G {0,1,2,.. .}, with u(N) + s(N) = n, 

• a homeomorphism c N : R™ — > R" = x W {N \ such that 



1 


1 








1 


1 


1 








1 


1 


1 



A = 








































c N (\N\) = B U (N)(Q,1) xB s (N)(0,l). 
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We denote such a quadruple by N. We usually drop the bars on the support and 
refer to it as N. Furthermore, 

JV C = 5^(0,1) x5^(fl,l), 
N~ = dB^(0, 1) x B^(0, 1), 
= B^)(°> x ) x 9 ^^(0, 1), 
N~ = cJ f 1 (N-), N+ = c N \N+). 

u(N) and s(N) are the nominally unstable and stable directions, respectively. 
The idea of a covering relation between two h-sets is that the image of the first should 
be mapped transversally across N~ (the exit set of N) and inside of N + (the entrance 
set of N). We formally define it. 

Definition 6.2 Assume N, M are h-sets, such that u(N) = u(M) = u and s(N) = 
s(M) = s. Let f : N — » R n be a continuous map. Let f c = cmo/oc^ 1 : N c — > R u x R s . 
We say that 

N M 

(N f— covers M) iff the following conditions are satisfied 

1. There exists a continuous homotopy h : [0, 1] x N c — > R u x R s , such that the 
following conditions hold true 

ha — fc, 

&([o,i],JV-)nM c = 0, 

h([0,l],N c )nM+ = fr 

2. There exists a linear map A : E" — > R u , such that 

h 1 (p,q) = (Ap,Q) for p£ S^(0, 1) and q e B~ s (0, 1), 
A(dB u (0,l)) Cl^^O.l). 
hi is called a model map for the relation N M. 

The maps that we study in the present paper are reversible, which we use to 
reduce the amount of computation. For such maps the following definition is useful, 
see also the discussion in (Zgliczyhski and Gidea 2004). 

Definition 6.3 Let N be an h-set. We define the h-set N T as follows: 

• The compact subset of the quadruple N T is the compact subset of the quadruple 
N, also denoted by N, 

• u(N T ) = s(N), s(N T )=u(N). 

• The homeomorphism c n t : R n E™ = R u( - nT ^ x M s ( ArT ) is defined by 

c n t(x) = j(c N (x)), 

where j : E"W x R s( - N ^ -» E S ( JV ) x R U W is given by j(p,q) = (q,p). 

Definition 6.4 Assume N,M are h-sets, such that u(N) — u(M) = u and s(N) = 
s(M) = s. Let g : fl C E" -» E n . Assume that g^ 1 : M — > E™ is well defined and 

continuous. We say that N d= M (N g—backcovers M) iff M T ^=4> N T . 
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Definition 6.5 Let N and M be h-sets. We say that N generically f -covers M 
(N 4=> M) ifN=^M or N 4= M. 

The main property of covering relations is contained in the following theorem 
(Zgliczyhski and Gidea 2004, Corollary 7). 

Theorem 6.6 Assume that we have the following chain of covering relations: 

N Q 4^ Nj. 4^ N 2 <=> . . . 4^ N k , 

then there exists a point x £ intATo, such that 

fi o fi-i o ... o f t (x) G intiVj, i = 1,. .. , k. 

Moreover, if — N , then x can be chosen so that 

fk o f k -! o ...o = x. 

In our proofs the hypothesis of Theorem |6.6| is verified using the routine 
checkCoveringRelations of (Progsl 2009). 

6.2. Cone conditions for h-sets 

Thcorcm |6.6| gives a computational tool to prove the existence of orbits with prescribed 
symbolic dynamics. To prove that such orbits are unique one would ideally require 
hyperbolicity of the map in a neighbourhood of the orbit. Typically, this is proved 
by constructing invariant cone fields. An alternative method to prove uniqueness is 
provided by covering relations with cone conditions, first described in (Kokubu et 
al 2007), the method is studied in further detail in (Zgliczyhski 2009), which we follow 
below. 

Definition 6.7 Let N C W 1 be an h-set and Q : R™ — > K™ be a quadratic form 

Q((x,y))=a(x)-(3(y), (x, y) e M" x M s , (34) 
where a : R"^) — > R, and (3 : — > R are positive definite quadratic forms. 
The pair (TV, Q) is called an h-set with cones. 

Definition 6.8 Assume that (JV, Qjv) and (M,Qm) are h-sets with cones, such that 
u(N) = u(M) = u, and let f : N -> R dm ( M ) be continuous. Assume that N M. 
We say that f satisfies the cone condition (with respect to the pair (N,M)) iff for any 
Pi,P2 € N c , pi ^ p 2 holds 

QmUM - fM) > Qn(pi - p 2 ). (35) 

Definition 6.9 Assume that (N, Qjv) an d (M, Qm) are h-sets with cones, such that 
u{N) = u(M) = u and s(N) = s(M) = s, and let f : N -> R M+S be continuous. 

Assume that N <= M. We say that f satisfies the cone condition (with respect to the 
pair (N, M) ) iff for any q\,qi € M c , qi ^ g 2 holds 

Qm(qi - q 2 ) > QN{f-\qi) - fc 1 ^))- (36) 
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Remark 6.10 Note, for details see the discussion in (Zgliczynski 2009), that stable 
manifolds are propagated backwards through a chain of covering relations with cone 
conditions. Similarly, unstable manifolds are propagated forwards through a chain 
of covering relations. In the planar case, used in the present paper, the Lipschitz 

constant of the stable manifold is given by \J~^, where alpha and beta are the constant 
coefficients of the corresponding quadratic forms. 

To verify that the cone conditions hold, we use the following Lemma (Zgliczynski 
2009, Lemma 8). 

Lemma 6.11 Assume that for any B G Dt c (Nc), the quadratic form 

V{x) = Qm(Bx) - Q N {x) (37) 
is positive definite, then for any p\,p2 G N c such that pi ^ P2 

Qu(fcipi) - /c(pa)) > Qn(pi - pa). (38) 
In our proofs the hypothesis of Lemma |6.11| is verified using the routine 
checkConeConditions of (Progsl 2009). 

Remark 6.12 The proof of Lemma \6 . 1 1\ also demonstrates that 

Qm(/c(Pi) " /c(P2)) - Qn{pi - Pa) = 0(\pi ~p 2 \ 2 ), 

and yields a uniform lower bound on 0(\p\ — P2 1 2 ) for each pair of h- sets: 
Qm(Mpi) ~ /c(Pa)) - Qn(pi -Pa) > £mn\pi ~P2\ 2 - 

Theorem 6.13 Let A = (£iij)i,7=o be a k x k matrix whose entries are either or 1. 
Assume that the following set of covering relations with cone conditions holds 

Ni<^>Nj ifdy=l. 



Then, each sequence in {0, 1, . . . , k — 1}^ (see Definition 3.4) is realised by a unique 
orbit. 

Proof. Assume that p\ S N Uo and p2 G A^ D , are such that p\ ^ P2, f l (pi) £ N Ui and 
f l {p2) € N iAji , for all i G Z, where ui — (..., U-i, (Jo, <^i> — ) € {0, 1, . . . , k — Let 



e = min^jj-w w =u tNjNn where etfjNi is as in Remark 6.12 There are three cases. 

(i) Let Qn^, (pi — Pa) > 0, then we have that for any i > 

Qn^PM - /c(Pa)) > QN^if^iPi) ~ /r X (P2)) + e\ft X (Pi) ~ /r X (Pa)| a 

> • • • 

> Qn uo (Pi - Pa) + e £ I fc (Pi )~fc(P2)\ 2 - (39) 

fc=0 

Since by assumption f l (pi) — P{P2) stays inside a compact neighbourhood of the 
origin, for all i G Z, where Qn^. is bounded from above, the left hand side is 
bounded. Hence fl(pi) — P C {P2) — ► 0, and therefore Qjv^. (/c(pi) — /c(Pa)) — * 0, but 
Qn^ (/c(Pi) - PApi)) > Qn„ q (Pi - P2) > 0, since Pi^p 2 - This is a contradiction. 

(ii) Let Q Na , (Pi - P2) = 0, then Qjv Wl (/ c (pi) - /c(Pa)) > 0, and the above 
argument holds for the pair / c (pi) and f c {p2)- 

(iii) Let Qn^ (pi ~ Pa) < we get the same contradiction for i < and ">" 
replaced by "<" . 



□ 
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Figure 4. The homoclinic tangle at the fixed point po- The stable manifold is 
given in blue, unstable — in red. The black boxes correspond to the /i-sets used 
in the proof. 



7. A transversal homoclinic orbit for F* 

The proofs of the theorems in this and the next sections were implemented in a 
computer program (Progsl 2009) using the software package CAPD (CAPD 2009) for 
interval arithmetic and covering relations. All computations were performed on an 
Intel Xeon 2.3GHz, 64bit processor with 16Gb of RAM. The program was compiled 
with gec, version 4.1.2. 

Numerical experiments indicate that for all F £ F* there exists a transversal 
homoclinic point p^, approximately located at (1.067707,0), to the fixed point po, 
approximately located at (0.577619,0), see Figure [3] We construct a sequence of 
covering relations and quadratic forms, and verify that the hypotheses of Theorem |6.6| 



and Lemma 6.11 are satisfied. Here, and in the rest of the paper, we use formula ( 11 1 
to provide bounds on the derivative of F. 

Theorem 7.1 Any map F € F* possesses a point p^ transversally homoclinic to p . 
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The map F 12 admits a locally maximal invariant set Cp\i 3 p*,f>rti.' 

Fl2 \c F i2 , ^2 |{0,1}-- 

homeo 

Proof. To avoid the "flip" of the system we study the second iterate of F, following 
the stable manifold from the homoclinic point to the fixed point. We construct 
four h-sets, Bi, i — 1,2,3,4, as in Table [T] box _B 4 around po, box B\ around 
Prh and two intermediate boxes B2, and B3. B4 is spanned by the unstable and 
stable eigenvectors of DF(po). The other boxes are experimentally rotated until two 
sides are roughly parallel with the stable manifold. We prove (using the routines 
checkCoveringRelations and checkConeConditions in (Progsl 2009)) that the 
following covering sequence with cone conditions holds 

Bi => B 2 => B 3 => B 4 => B^ 

The map F is reversible, by reversing time and reflecting the h-sets in {u = 0} our 
cover is also a back-cover from p^ to po along the unstable manifold of pq. Therefore 
we also have the following covering relations with cone conditions 



Bi T(B 3 ) T(B 2 ) B t . 



Existence of a homoclinic orbit follows from Theorem 6.13 The transversality of the 



intersection follows from the fact that we use cones with Lipschitz constant less than 
one, Lip = 0.95. 



□ 

Note that both po and p^ depend analytically on F. Part 1 of the first Main 
Theorem follows immediately, since F* £ F t . 



Box # 


Centre 


"Unstable" Direction, e" 


Bi 


(1.067744,0.0) 


(0.271460881389316,0.962449473933741) 


B 2 


(0.90733000, -0.328353) 


(0.653774812595729,0.691649449891044) 


Bi 


(0.66401000, -0.0727302) 


(-0.736700425093906, -0.676219257095255) 


B 4 


(0.5776305,0.0) 


(0.788578889012330,0.614933602760558) 



Table 1. Boxes Bi used to prove the existence of a transversal homoclinic orbit. 
The size of each side of B4 is 0.035|e4|; the size of Bi, i = 1,2,3, is 0.01|e u |. 
The "stable" spanning vector is defined as e s = T(e u ). The covering sequence is 
B 1 => B 2 => B 3 => Bi B 4 . 



Remark 7.2 The numbers describing boxes B 4 in the Tableware non-representable 
on the computer. We have used these numbers as the input for the routine 
checkConeConditions. During its execution, the numbers are rounded up to the 
nearest representable number; therefore, strictly speaking, Theorem 7.1 has been proved 
for boxes different from (although very close to) those reported in the table. Since the 
specific values of the true parameters are largely irrelevant, we will not go into pains of 
computing and reporting bounds on these parameters. This will be the situation through 
out the paper: all specific numerical data appearing in various tables and statements 
of theorems are (very accurate) approximations of the true data for which the results 
have been proved on the computer. 
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8. A bi-infinite heteroclinic tangle for 

In this Section we will demonstrate that the fixed point has a sequence of hyperbolic 
periodic points {z k }, fceZ, 



and 



k — > — c 



such that their stable and unstable manifolds form a heteroclinic tangle. 

We first demonstrate that there exists a heteroclinic orbit between the fixed point 
zq = po and the period two point z\ = A»(po)- 

We proceed as in the previous section: we construct a covering sequence of /i-sets, 
as in Table [2] consisting of: one box Bq around po, one box B\ around A*(po) and 11 
intermediate boxes on the orbit. In this sequence, Bq is defined as A~ 1 (Bi). 



Box # 



Centre 



If 



I" 



Bi 
B 2 
B 3 

Br, 

B 6 
B T 
B 8 
B 9 
Bio 



(-0.143816,0) 0.00072 0.00072 

(-0.141389,-0.000460854) 0.00054 0.00015 

(-0.133448,-0.00190821) 0.00045 0.00015 

(-0.098565,-0.00713661) 0.00047 0.0002 

(0.0582909, -0.00620757) 0.00055 0.00015 

(0.106679,0.0110872) 0.0004 0.0001 

(-0.2842,0.070894) 0.0004 0.0001 

(0.533255,0.0331362) 0.0003 0.0002 

(0.575172,0.00190701) 0.0003 0.0002 

(0.577516,0.000132348) 0.0005 0.0002 

(0.577619,0) 0.001 0.001 



(0.982114,0.188285) 
(0.982114,0.188285) 
(0.982114,0.188285) 
(0.987,0.16072) 
(0.994,-0.10938) 
(0.975,-0.222205) 
(0.84,0.542586) 
(0.8, -0.6) 
(0.788579, -0.614934) 
(0.788579, -0.614934) 
(0.788579, -0.614934) 



(0.982114,-0.188285) 
(0.982114,-0.188285) 
(0.982114,-0.188285) 
(0.997, -0.0774016) 
(0.995, 0.0998749) 
(0.982114,0.188285) 

(0.866025,-0.5) 
(0.788579,0.614934) 
(0.788579,0.614934) 
(0.788579, 0.614934) 
(0.788579,0.614934) 



Table 2. The data used to prove the existence of a heteroclinic orbit for F% 
between zq and z\. Vectors e"' s are the "stable" and "unstable" spanning vectors 
of the rectangles. The length of the sides of the rectangles Bi is 2 ■ i"' s • |e"' s |. 
The covering sequence is B\ => B\ => B2 => B3 =>■ B4 => B$ =>■ Bq => Br => 
B 8 => B 9 => B 10 => B11 => B Q = A.^Bi) =4> B = A* ^Bi). The first and the 
last two covering relations are for F%. the other ones — for Fr. 



Since the orbit is on the tangle between the unstable and stable manifolds, 
there are parts of it where it is hard to construct expanding and contracting 
directions directly for the second iterate of the map. The fourth iterate, however, 
exhibits nice hyperbolic like behaviour along a pseudoorbit. We use the programs 
checkCoveringRelations and checkConeConditions to verify that the hypotheses 



of Theorem 6.6 and Lemma 6.11 are satisfied for on the sequence of /i-sets. In 
addition, we prove that they are satisfied for F% for the covering B\ =>• B\ (which 
implies that the covering relations and cone conditions are also satisfied for Bq =>■ Bo). 
We use cones with Lipschitz constant 0.9 in the proof. 

The reversibility of F„ implies that the existence of a heteroclinic orbit from Z\ to 
zq entails the existence of a heteroclinic orbit from zq to z\ . It also implies existence 
of homoclinic orbits for both points, that enters a small neighbourhood of the other 
one. The sizes of these neighbourhoods are given by the sizes of the boxes Bq and B\ 
in Table 

To construct the pseudoorbit we non-uniformly distribute points in the unstable 
eigendirection of z\, and locate an initial point that passes as close as possible to z\ 
and zq. The finest discretization of the approximate unstable vector for the period 
two point is done on the scale 10 -14 . The first point on the pseudoorbit is in the 
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unstable eigendirection of z\ at a distance of 0.00058120854283, and the last point is 
at a distance of 0.000167835 from Zq. 

We will now demonstrate the following eight covering relations: 

Bi => B , Bi B\\ B a <= Bi, Bi <= B\\ 

Bo =^> Bo, Bo <= Bo; Bo ==> B\, B\ <= Bo- 

Step 1. Set 

iB? = F~ 4 (F 4 (Bi) n F~ 4 ( 

. . . n F" 4 (F 4 (B 8 ) n F" 4 (F 4 (B 9 ) n F" 4 (F 4 (B 10 ) n B n ))) . . .)), (40) 
1 b\ = F- 4 (i?(Bi)nF- 4 ( 

. . . n F-^FtiB!) n f- 4 (f 4 (B0 n f-^f 4 ^) n Si))) . . .)) (41) 

Clearly, the iterate F* = F 42 is well defined and analytic on iB° U iBj, and 

iB? ^ B , iBi 1 ^ Bl (42) 

The set B\ is a parallelogram. Let 7rf be a projection on the span of ef "along" 
the vector e", and 7r" be a projection on the span of e" "along" the vector e\ . Then 
the sets \B\, i — 0,1, are "vertical" in the sense that nl(iB^) — nl(Bi), while the 
sets F t (iB{), i = 0, 1, are "horizontal" in the sense that ttqF^^B^) = ttq(B ) and 
ir?F.( 1 Bl)=n?(B 1 ). 

Step 2. Since the set B\ is symmetric: T(B\) = B\, we have that 

T( 1 B 1 °)^B , TdB^^Bx. (43) 
Step 3. We get from the period doubling equation: 

F,(A; 1 ( 1 B 1 1 ))^A; 1 (F 2 ( 1 B 1 1 )). 

Clearly, since F* is well-defined and analytic on iB\, and \B\ ^> \B\ there exists a 
vertical subset 2 B\ C \B\ on which F 2 is analytic. Set iB§ = A^ 1 ( 2 Bj). Then, since 
Bo = A-^Bi), iB ° covers B : 

xBg^Bo, T(iB°)^fl . (44) 

Siep ^. Let iBj D iB^ be any set such that e < dist^^^ofe iB") < 2e for 
some sufficiently small e > 0, chosen so that F* is analytic on iBj. The set 
xB 1 = F" 1 (TdB?)) r|B satisfies 

iBo 1 ^> Bi, F(iBq) ^ Bx. (45) 

Set 

Cf = lim C| , 

- r * fc— »oo ^* 

where 

c*. = f^c*; 1 ) n f- 1 ^; 1 ), = K(ul =0l B{) n f- 1 (u 1 1 j=0 t( 1 b|)). 

We have proved the following: 
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Lemma 8.1 The map = Ff 2 admits a locally maximal invariant setCp on which 
its action is homeomorphic to the full Bernoulli shift on {0, 1} Z . 

Proof of Main Theorem [7] ( 2-4): 

The covering relations (42 1 imply that, for any to, there exist vertical sets 

m Bi C 1B1 and m B\ C \B\ such that F™ is analytic on them, and m B\ ==>■ Bq 

su\& m B\^B x . 

Consider the map F 2 for some integer k > 1. The fixed point equation 

A~ fe o Pf o = I", (46) 

implies that 

K\T{ n2k Bl)) % K\T{B )) = A~ k ~ 1 (T(B 1 )), 

K k {n2«B\) S A,- fe (Bi), 

K k {T{ n2k B\)) S A-'iB,). 

Below we will use the notation that A =^> B if there exist some vertical set 
A C A (depending on n and A) on which F™ is analytic, and A B. A ■(= B and 

A <J=S» B are used similarly The covering relations ( 42 ) together with the fixed point 
equation (46 1 imply that, in this notation, we get for any integer k > 

A k (B x ) A k -\B t ) . . . ^ By, 

Aj(T(Bi)) <^> t£-\T{B{)) K <U T(B 1 ). 

We have, therefore, shown that for any two natural N and M, there exists an 
integer n = 2 N , such that 

Af (B 1 ) S> A?~\ Bl ) &...& Bl & A-\B,) M> 

A~ M +\ Bl ) & A~ M ( Bl ) S> A~ M +\ Bl ) M> 



=£» A-^B,) ^B^...^ Af- 1 ^) 4^ Af 



Set 



Aat,m — [J A* (0( fc )S°) U Aj; (^Bl) , 

fe=-Af 

where 0(fc) = 1 if k > and 0(fc) = 2 Ar_fe if k < 0, and define recursively: 



4„ ee ^"(A^m) n F- n (T(A NM )) and Z)„ ee F^Z'T 1 ) n F^Z*" 1 ). 

Clearly, the set Zj™ ee Z|?„ is a locally maximal invariant set for F™ , appearing 
in part 3) of the Main Theorem 1. 
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F™\z F n is homeomorphic to a topological Markov chain, which is mixing. This 

together with dist(A~ fc (i?i), 0) 00 an d dist(A»(Bi), 0) fej= ^ implies the 

existence of unbounded and oscillating orbits. This proves part 4) the Main Theorem 
1. 

Finally, we notice that if Oq is a heteroclinic orbit between points z\ and zq for 



the iterate F* = (or vice versa), then the fixed point equation (46 1 implies that, on 
one hand O k = A* (Co) is a heteroclinic orbit for iterate between periodic points 
z k+ i and z k : 

z k = A k (z ); 

on the other hand if p G 0- k = A^ fc (O ) then A*(p) G O , and so is q = F^*(A*(p)), 
and 

F.(p) = h- k {F? (A*(p))) = A- k (q) G 0_ fe . 

Therefore, contains a heteroclinic orbit between fixed points z%- k and 
for iterate F*. Stable and unstable manifolds of points 2^ and Zfc+i intersect in a 
two-sided tangle, see Figure [5] This proves part 2) of the Main Theorem 1. 

□ 

9. Distortion tools 

We will say that a pair (/C, ip) is a dynamically defined Cantor set if /C C i? is a Cantor 
set and -0 : /C 1— > /C is a locally Lipschitz expanding map, topologically conjugate to 
the full Bernoulli shift a : {0, 1} N h-> {0, 1} N . 

A Markov partition V = {/Coj^-i} of {K,,ip) is a partition of /C in two disjoint 
subsets, K. = /Co U fCx, /Co fl /Cj = 0, such that t/^/q : /C; i— > /C is a strictly monotone 
Lipschitz expanding homeomorphism. 

Given a sequence (a , . . . , a n _i) G {0, 1}™, denote 

n-l 

/C(a ,...,a n _i) = P| tp~ l (K, at ). 

i=0 

A bounded component of K \ /C is called a gop of JC. Gaps can be ordered as 
follows. The unique gap of order is the interval Uq = K \ (jCo U ICij , where A 
denotes the convex hull of a set A C M. A connected component of 

K.\ [J JC(a , . . . , a n -i) 

(o ,-,ai)e{0,l}« 

is called a gap of order n if it is not a gap of order k < n. Any gap of JC is a gap of 
some finite order. 

Definition 9.1 Given a gap U of K, of order n, we denote Ljj (respectively, Rjj ) the 

unique left (respectively, right) adjacent to U interval of the form /C(ao, . . . , a„_x), 
(a ,...,a n _!) G{0,1}". 
The numbers 

t l (K) = inf I M : U is a gap of £ j , (47) 
t r (K) = inf I l^i : U is a gap of /c| (48) 
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T A, 



b) 




c) 




Figure 5. Intersections of a) the stable manifolds of z\ and F*(zi) (blue) with 
the unstable ones for z%, F,(z-2), F%(z2) and Ff(z2) (red and magenta); b) the 
stable manifold of zg (blue) with the unstable ones for z\ and F t (zi) (red); c) 
the stable manifold of Z—\ (blue) with the unstable one for zq (red). 
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will be called the left and right thicknesses of K. The numbers 

r L {V)^^4 and r R {V) ee ^ 



\U \ ^ ' ~ \U \ 

are called the left and right thicknesses of the Markov partition V . Given a Lipschitz 
expanding map i/i:/nl,/cR, define the distortion of on I: 

( \qh(y) - lh(x)\\z - Zll 

Distty, I) = sup { ^, z£/: ^ x ^ x} log | _ ^ y _ g | | e [0,oo] (49) 
( notice, ip(z) ^ ip(x), ip(y) ^ tp( x ) )■ 

Definition 9.2 The distortion of the dynamically defined Cantor set (JC,t/j) is 

Dist^/C) = sup (o0i ..., 0n _ 1 ) £ {o,i}" Di st(V' n ; ^(a , . . . , a n _i)). (50) 

Lemma 9.3 (see fPa/is and Ta/cens (Duarte 2000)) Let (/C,V) &e a 

dynamically defined Cantor set with a Markov partition V and distortion Dist^,(/C) = c. 
Tften 

e- c r L (P) < t l {K) < e c r L (V), e~ c T R {V) < t r {K) < e c r R (V). (51) 

This lemma is important since it allows one to estimate the thicknesses of a Cantor 
set in terms of thicknesses of its Markov partition. 

Definition 9.4 (see (Duarte 2000), (Gorodetsky and Kaloshin 2008)) Define T to be 
the set of all maps F : A U Ai R 2 such that: 

1) A and Ai are compact sets, diffeomorphic to rectangles, with non-empty interior; 

2) F is C 2 on a neighbourhood of A U Ai, that maps A U Ai diffeomorphically 
onto its image; 

3) the locally maximal invariant setCp = n„ e z-F 1 ~"(AoUAi) is a hyperbolic set, and 
the action of F on Cf is conjugated to the Bernoulli shift a : {0, 1} Z {0, 1} Z . 

4) V = {A , Ai} is a Markov partition for F : Cf i— > Cf, in particular, F has two 
fixed points po G A and p\ e Ai, whose stable and unstable manifolds contain 
the boundaries of A and Ai. 

We will now introduce our main tool for computing the Hausdorff dimension 
the Duarte Distortion Theorem (see (Duarte 2000), (Gorodetsky and Kaloshin 2008)). 

Definition 9.5 Given positive constants C, e and 7 define T(C, e, 7) to be the class 
of maps F : A U Ai ^ R 2 , F E T, such that 

1) diam(A U Ai) < 1, diam(F(A ) U F(Ai)) < 1; 



2) the derivative of F, DF(x,u) 
satisfies on Aq U Ai 

a) detDF(x,u) = I, 

b) \d\<l< \a\ < f , 

c) |&|,|c|<e(|a|-l); 



a(x,u) b(x,u) 
c(x,u) d(x,u) 



, where a, b, c and d are C 1 , 
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3) functions a = aoF~ l , b = bo F~ l , c = coF^ 1 , c = co F^ 1 on F(A ) U F(Ai), 



appearing in DF (x,u) 



satisfy 



a) 


db 




dd 




db 




.95 




dd 




dc 


dx 




du 
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du 


j 


dx 
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dx 




dy 


b) 


da 




db 




db 




dc 




dc 




dd 


du 




dx 
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du 
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dx 
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du 




dx 


c) 


da 
du 




dd 

dx 


<7|5|( 


|5|- 


-1), 






d) 


da 
dx 




dd 
du 


<7|5|( 


|a|- 


-i); 







< 7(|a| -(m 

<7(I«K55) 
(54) 
(55) 

4) the variation of log \a(x, u)\ in each rectangle Aj is less or equal to 7(1 — oti), 
on = max( X)U ) SAi |o(a;,u)|; 

5) the gap sizes satisfy: 

dist(A ,A 1 ) > -, dist(F(Ao),F(A 1 )) > -. 

7 7 

Remark 9.6 Unlike in (Duarte 2000) and (Gorodetsky and Kaloshin 2008) we do 
not require the eigenvalues of the maps F G T at the fixed points po and p\ to be 
positive. The positivity of the eigenvalues is important only in the construction of the 
boundaries of the two components A and Aj of the Markov partition V (see Lemma 



10. 5| for the construction of V in our case), but does not enter the proof of Theorem 



We will now continue with the notion of stable and unstable Cantor sets. Denote 
the stable and unstable foliations of Cf as J- s and T u : 



T s = {connected comp. of W S (C F ) n (A' U A^)}, 

T u = {connected comp. of W U {C F ) n (F(A' ) U F(A[))}. 



(56) 
(57) 



Also, define 



I s = W£, c (po) n A„ and l u ee WZJfr) n A . 
Then, the Cantor sets 

JC S ee C F n I u and JC U = C F n I s 



(58) 



can be identified with the set of stable leaves of J- s , respectively unstable leaves of 
T u . Define the projections 7r s : Cf >— ► IC S and ir u : Cf K. u in the following way: if 
p € Cf then 

7r s (p)=Wf oc (p)nT" and tt u (p) - W£ c (p) n I s , 
and the maps ip s : K. s ^ K. s and ^" : /C" 1— > /C" as 

i/> s ee tt s o F and ^l" = 7T U o F~^ . 

The pairs (IC s ,ip s ) and (/C u , V 1 ") are dynamically defined Cantor sets with Markov 
partitions 



V s ee {I u n Aq,2 u n Ai} and P" ee {X s n F(A ),X S n F(Ai)}. 
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Theorem 9.7 (Duarte Distortion Theorem, (Duarte 2000)) If F £ T{C, e,j), then 
the distortion of the dynamically defined Cantor sets (K. u ' s , ip u ,s), i>u,s — t«, s ° F, is 
bounded: 

Disfy Ui . (JC U ' S ) < D{C, e, 7 ) = 4(C + 3) 7 + 2e. (59) 

in particular, 

e -D(C,e,7) Ti (7?«.«) < t l (/C u ' s ) < e 2, ( c,£>7 Vi,(7" 1,a ), 
e- D ( c ' £ >T , )T B (7^ ,a ) < t^(/C"' s ) < e D ( c ' £ ^r fi (P tt ' s ). 

The Duarte Distortion Theorem has been used as a powerful tool in conservative 
dynamics on several instances. It has been proved, and applied in (Duarte 2000) 
to demonstrate that accumulation of a locally maximal invariant set by periodic 
elliptic points is generic in two parameter families. Furthermore, (Gorodctsky and 
Kaloshin 2008) uses it to show that generic unfoldings of homoclinic tangcncies of two- 
dimensional area-preserving diffcomorphisms give rise to hyperbolic sets of dimension 
arbitrarily close to 2. To show this, the authors of (Gorodetsky and Kaloshin 2008) 
complement the Duarte Distortion theorem with the following results, which will be 
useful to us: 

Lemma 9.8 The Hausdorff dimension of a Cantor set JC satisfies 

dim#/C > d, 

where d is the solution of 

r L (IC) d + T R (JC) d = (1 + tl{K) + T R (JC)) d . 



Lemma 9.9 Suppose that for given th > 0, t R > the solution d of the equation 

t d L + t d R = (l + t L + t R ) d 

is in (0, 1), then 

'log(l + T ^-) log(l + T ^- 
d>max<; ) ) ^fV. (60) 



log (l + i±^) ' log (l + i±^) J 



10. A horseshoe for 

In this Section we demonstrate the existence of a horseshoe for the third iterate of 
GeG,. We start with definitions. 

Definition 10.1 (Full component) Suppose AcMcM" = E^elR' is homeomorphic 
to a rectangle (a set of the form D\ x D 2 C R k © M 1 ), and let F : U ^ R n be a 
diffeomorphism. A connected component A = F(A ) o/Afl F(A) is called full, if 

1) V 2 (A ) = D 2 , 

2) for any z G A , "Pi| F (A n(D 1 x'P 2 ( z ))) is a bijection onto D 1 . 

Here, V\, 2 = V\^ph, h the homeomorphism A ^ DixD 2 , V\ andV 2 are the canonical 
projections on R k and R l . 
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Definition 10.2 (Two-component horseshoe) Let U C K™ be an open set, then the 
set A C hi homeomorphic to a rectangle D\ x D 2 is called a two- component horseshoe 
for the diffeomorphism F : U i— > M. n if A n F(A) contains at least two full components 
Ao and A\ such that, 

1) "P 2 (A U Ai) c intL> 2 , Vi(F^ 1 (A U Ai)) c intD i; 

2) -D-F 1 |f- 1 (A uAi) preserves and expands an unstable ("horizontal") cone family on 
F-^AqUAi); 

2) -Di 7 ^ 1 |a uAi preserves and expands a stable ("vertical") cone family on AoUAi. 

By a standard result of the theory of dynamical systems, if A is a two-component 
horseshoe for F, then the action of F on the locally maximal hyperbolic set Cp = 
n^_ oc F l (A) is homeomorphic to the Bernoulli shift a 2 on {0, 1} Z . 

To demonstrate the existence of a horseshoe for the third iterate of G G G,, we 
construct two rectangles: A' which contains the fixed point po and A^ which contains 
the period three point pi. We will prove that these two rectangles constitute two full 
components of a horseshoe, and are such that G(A' Q ) n Wf oc (j?o) and G(A' 1 ) nWf oc (po) 
contain the components of the Markov partition of the locally maximal hyperbolic set 

£?■ 

The centres of the rectangles and sizes of the spanning vectors, are given in 
Table [3] Note that these boxes are symmetric with respect to the involution 
T : (x,u) i ► (x, — u). Therefore, we only have to prove that the covering relations 
hold for G*. 

Theorem 10.3 For any G G G*, the inverse map G _1 admits a horseshoe whose 
two full components are the sets Aq and A' l7 defined in Table ^ A' Q and Aj are 
independent of G. 

Proof. To demonstrate the existence of the horseshoe we first use the routine 
checkCoveringRelations to show that 

A' A' iS> A' =^5> A' A' 



where the /i-sets are as in Table [3] This implies that 1) in the Definition 10.2 holds. 

The remainder of the proof of the existence of a horseshoe is done for the map 
G- 1 = T o G* 1 o T~\ where T(x,u) is some coordinate transformation that 
approximately diagonalizes the derivative DG 4 on A' U A^. (See (Progs2 2009) for 



Component 


Centre 


"Stable" Scale 


"Unstable" Scale 




(0.670198,0.0) 


0.083 


0.083 


a; 


(-0.441811,0.0) 


0.0655 


0.0655 



Table 3. The rectangles that approximate the Markov partition for 
the horseshoe of G*. The spanning vectors are approximate stable and 
unstable vectors for DG*(po,i). The "stable" vectors for po and pi, 
respectively, are eg = (0.788578889012330,-0.614933602760558) and ef = 
(0.750925931392967773,0.660386436536671957). The "unstable" vectors arc 
e = ^( e o) an( ^ e l = ^X e i)- The length of the sides of the rectangles Aq 
and A' x is 2 ■ stable/unstable scale • |gq | . The covering sequence, which holds 
for G*, is A?, =>• A;, => A' => A' AL 
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programs used in this part of the proof). We have used the following coordinate 
change: 

T(x, u) = 0.55 (x - 0.5, u + 0.0265461632116977382 

- 0.658388496175704694(x - 0.577619) 

- 0.611583723237529069(:r - 0.577619) 2 

+ 0.102408008658008658(a; - 0.577619) 3 ). (61) 

Denote Aq 1 = T(A' Q t ). We show that there exist an invariant unstable cone 
field C u on A' U A' x , 

C u (p) = {v=( Vl ,w 2 )e]R 2 :-0.15<^<0.2}, p £ Aq, (62) 
C;(p) = {v=(t. 1 ,« 2 )eR 2 :-0.1<-<0.05}, p £ A[. (63) 

Vl 

We have verified that 

DGt,(p) ■ v G int C„(G* (p) ) , for all v G C u (p) , p £ A' Q U A[ , (64) 
and that vectors inside the cone field are expanded: 

\DG*(p) ■ v| = A|v|, v e mtC u (G4p)), 

where the expansion rate A is in A = [A_, A + ] 

A^ = 6.091, A+ = 26.892. (65) 

Verification of invariance and expansion of the cone fields is carried out (among 
other things) in the routine Distortion of (Progs2 2009). 

Reversibility and symmetry of the partition imply that there exists a stable cone 

field 

c s {p) = i?r(T- 1 (p))(r( J Dr- 1 (p)c„(p))) (66) 

on (^(Aq U A[) n ^Ag U a0, which is mapped into its interior and expanded with 
the rate A by DG 



-l 



□ 

Lemma 10.4 The cone fields 

C u ( P ) = DT- 1 (T(p))C u (T(p)), peA^UAi, (67) 

C s (q) =T(C u (q)), q £ G*(Aq U A^) n (A' U A^) , (68) 

are transversal in the sense that the angle between any u G C u (q) and v G C s (q), 
q £ G*(A' U Ai) n (A' U A[) is bounded from below. 



Furthermore, the leaves of the foliations T s and T u , defined in (56 1 — (|57|, are 
graphs over the x-axis. 

Proof. We verify on the computer (this is done in the routine Distortion in 
(Progs2 2009)) that any v G L»r- 1 (T(p))C r u (r(p)), p £ A' , satisfies: 

— > 0.461, 

Vl 
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and that any v e DT- 1 (T(p))C u (T(p)), p € A' x , satisfies: 

— < -0.675. 

Vi 

This implies that the angle between any u G C u (q) and v £ C s (q), q € G*(A U A[) fl 
(Aq U A' x ) is bounded from below. 

□ 

The horseshoe for G is illustrated in Figure [6j for both original (a) and T (b) 
coordinates. 

Note that the full components Aq, and A^ are chosen uniformly for all G € G*. 
This is not compatible with the definition of full components in Definition [9^4] We can, 
however, prove that any G € G* has full components, as in Definition |9.4| contained 
in the components constructed in Theorem 1 10. 3 1 

Lemma 10.5 For any G € G* there exist full components Aq C A and A\ C A x as 
in Definition\9.4\ 



Proof, (see (Progs2 2009) for programs) Let Ge G,, by Theorem |1 0.3 

Pi = T{Px) <s A' x n G(A'), 

hence there exists e > such that B e (px) C A'. Therefore, VV £ " = VV£ c (pi) ri-B e (pi) C 
A . Let t h-> VV"(t) be some parametrization of VV". For sufficiently small e the 
tangent vector r(i) to VV"(t) is close to e"(pi), that is, there exists e > such that 
r(t)€C u (W«(i)). _ 

The cone conditions imply that W" (t) is expanded in the forward (unstable) cone 
into cvxG(A ): there exists a positive integer N such that Gr (VV") D (A ^ ^ 
and G^CVV") D (A' 10 )2 ^ where (A 01 )~, i = 1,2, denotes the two components of 
the exit set of A x (see Definition |6.lj) . 

Set b Q ee G Ar (>V^)nA and 6i ee G A, (W e ")nA' 1 . G(6 ), stretches across A', by the 
set of covering relations in Table 3 1 that is G(6o)n(A 01 ) 7 " and G(6o)n(A 10 ) r " 7^ 0. 

We set t ee G(6 ) H A and i x ee G(6 ) H A' r 

Via a similar construction for G _1 , we can show that there is a pair of stable 
leaves h and n such that Zi n (A'^ ^ 0, i = 1,2, and n n ^ 0, i = 1,2. 

Here (A i = 1,2, denote the two components of the entrance set of A j (see 



Definition 6.1 1. Furthermore, there is a pair Iq and ro with similar properties. 

Since bi and t\ stretch across A x "horizontally", and li and r% stretch across A x 
"vertically" , bi intersects l\ and r\ at least once, and t\ intersects l\ and rj at least 
once. Clearly, p\ is one of these intersection points. Therefore, 61, t\, r\ and l\ bound 
an open set U £ A x with at least one connected component. We define Ai as the 
closure of the connected component that contains p± in its boundary. 

Ao is constructed in a similar way: it is the closure of the connected component 
bounded by bo, t , l Q and ro that contains po in its interior. 

□ 
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Figure 6. a) The two components Aq and A'j of the horseshoe map across 
A U Aj under G and G~ 1 . A Q is the diamond on the right, Aq — the diamond 
on the left. The entrance sets are given in red and magenta, the exit sets are in 
green and blue. The two sets whose longer sides in red and magenta are G(A' i ), 
i = 0, 1, the sets with longer blue and green sides are G~ 1 (A' i ), i = 0, 1. b) The 
same horseshoe in the coordinate T. The local stable manifold is given in cyan, 
the boxes signify the endpoints of the Markov partition for /C". 



11. Hausdorff dimension of the locally maximal invariant set for 



According to Lemma |9.8| and the Duarte Distortion Theorem |9.7[ the knowledge of 
the left and right thicknesses of the Markov partition for the unstable Cantor set K, u 
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together with the distortion of this Cantor set is sufficient to compute a lower bound 
on its Hausdorff dimension. 

Lemma 11.1 For any G G G* the intersections of the local stable manifold Wf oc {po) 
with 8G(Aq) and dG(Ai) are at the points fcg,/cQ and k\,k\, respectively. 
Where, 

k b G [0.041476215,0.04509142] x [0.01162711,0.015242315], 
fc* G [0.034578335,0.03698849] x [0.023247785,0.02565794], 

and 

k\ G [0.099814,0.1022241] x [-0.07865165,-0.07624155], 
k\ G [0.09487115,0.09660585] x [-0.06980765,-0.06710275], 

respectively. 

Proof, (see (Progsl 2009) for programs) We prove that there exists two points 
qi ,q 2 G W£ c (pi) n G(Ai) n W s ( Po ), such that 

G(<zi)nwf oc (p o )^0 

G(q 2 ) n A ,4 

G 2 ( (?2 )nA 1 ^0 
G 3 fe)nwf oc ( Po )^0 

We now describe the construction of gi. Numerical experiments indicate that there 
exists a heteroclinic orbit between p , and the period three point pi, approximately 
located at (—0.527156, 0). To prove this we construct a covering sequence of h-sets, as 
in Table [4] consisting of: one box around po, one box around p\ and five intermediate 
boxes on the orbit; two of them close to the unstable vector of p\ and three of them 
close to the stable vector of pq. We use the programs checkCoveringRelations and 
checkConeConditions to verify that the hypotheses of Theorem |6.6| and Lemma |6.11| 
are satisfied on the sequence of h-sets. 

We emphasize that the choice of boxes is very delicate. To construct the sequence 
we first construct relatively large boxes around pq and p\, spanned by their unstable 
and stable eigenvectors. Second, we reduce the size of these boxes until we are able 
to prove that the hypothesis of Lemma |6.11| is satisfied. Third, we non-uniformly 
distribute points in the unstable eigendirection inside of the box containing the period 
three point; on the smallest scale the points are only separated by 10~ n . Fourth, we 
iterate, non-rigorously, these points to locate a heteroclinic pseudoorbit, which starts 
within the box containing p\ and ends in the box containing pq. The first point on 
the orbit is in the unstable eigendirection of the period three point, at a distance of 
0.000158885 from our approximation of the period three point, and the last point is 
at a distance of 0.00177401 from the fixed point. Finally, we use the points on this 
pseudoorbit as centres of our covering sequence. The first two intermediate points are 
on the local unstable manifold of the period three point, and we use its unstable and 
stable eigenvectors to span the corresponding h-sets. The other three intermediate 
points are on the local stable manifold of the fixed point, and we use its unstable 
and stable eigenvectors to span the corresponding h-sets. Note that if we change the 
location of the original point on the stable vector by as little as 10 -8 , the proof fails. 
The existence of q\ follows from the set of covering relations constructed in Table [4] 
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To find <72 we do a similar non-rigorous computation to locate a pseudoorbit with 
the properties above. The corresponding sequence of h-sets and covering relations is 
given in Tabled Finally, we put k\ = TG(q\), fc* = TG 2 ( qi ), k b Q = TG 3 ( qi ), and k\ = 
TG 3 (q2). By the construction of the components of the Markov partition Ao and Ai, 
it is immediate that {k l , k b } = W{ oc (p ) n <9G(A ) and {k{ , k\} = VVf oc (p ) n 5G(Ai). 
Evaluating at q /i with G* yields the enclosures given in the statement of the 
lemma. 

□ 

Consider a G € G*, by Lemma |"l0.5| the two components of the Markov partition 
for K, u are contained in the intersections 

G(A ) n (Wf oc (p ) n A ) , G(A0 n (wf oc (p ) n A ) . 



Box # 


Centre 


Eigenvector 


Scale 





(0.577619,0) 


Fixed 


0.004 


1 


( -0.527156, 0) 


Period Three 


0.0005 


2 


(-0.524339,-0.00302553) 


Period Three 


0.001 


3 


(-0.458451,-0.0695428) 


Period Three 


0.002 


4 


(0.683671,-0.0907786) 


Fixed 


0.002 


5 


(0.565061,0.00974195) 


Fixed 


0.002 


6 


(0.578698,-0.0014086) 


Fixed 


0.003 



Table 4. The boxes used to prove the existence of a heteroclinic orbit between 
the period three point and the fixed point, passing through q\. Note that the 
numbers in this table are in the original coordinates. The point qi corresponds to 
box #3. The spanning vectors are the stable and unstable vectors for either the 
fixed point or the period three point. Note, the stable vector is constructed by 
reflecting the unstable vector in {u = 0}. The approximate unstable eigenvectors 
for the fixed and period three points are (0.788578889012330,0.614933602760558) 
and (0.682259082166558,-0.731107710380897), respectively. To construct the 
corresponding parallelogram the spanning vectors are scaled as indicated in the 
table. The covering sequence isl=>l=>2=>3=>4=>5=>6=>0=>0. 



Box # 


Centre 


Eigenvector 


Scale 





(0.577619,0) 


Fixed 


0.002 


1 


( -0.527156, 0) 


Period Three 


0.00015 


2 


(-0.527035,-0.000143356) 


Period Three 


0.0001(0,9) 


3 


(-0.523879,-0.00351678) 


Period Three 


0.0003 


4 


(-0.447416,-0.0799474) 


Period Three 


0.0003 


5 


(0.787595,-0.00888005) 


Fixed 


0.0005 


6 


(-0.36687,0.0109682) 


Period Three 


0.0003 


7 


(0.67407,-0.0819091) 


Fixed 


0.001 


8 


(0.566215,0.00882267) 


Fixed 


0.001 



Table 5. The boxes used to prove the existence of a heteroclinic orbit between the 
period three point and the fixed point, passing through qi- Note that the numbers 
in this table are in the original coordinates. k\ corresponds to box #7 and 52 to 
box #4. The spanning vectors are the stable and unstable vectors for either the 
fixed point or the period three point. Note, the stable vector is constructed by 
reflecting the unstable vector in {u = 0}. The approximate unstable eigenvectors 
for the fixed and period three points are (0.788578889012330,0.614933602760558) 
and (0.682259082166558,-0.731107710380897), respectively. To construct the 
corresponding parallelogram the spanning vectors are scaled as indicated in the 
table. The covering sequence isl=>l=>2=>3=>4=>5=>6=>7=>8=>0=>0. 
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Therefore, our next task will be two find an enclosure of the local stable manifold 
at the fixed point p a . 

By the property of h-sets with cone conditions mentioned in Remark |6.10| it 
is possible to propagate stable manifolds backwards through a sequence of covering 
relations. We use this property in the program encloseStableManif old to enclose 
the local stable manifold VV( oc (po) H Ao; estimates of the length of its intersection with 
G* (A) are needed to estimate the left and right thicknesses of the Markov partition 
for IC U , and ultimately the Hausdorff dimension. 

To enclose the stable manifold of the fixed point, we work in the coordinates (61 ). 

First, we choose 200 points q± distributed uniformly on the interval 

(p 0l p + 0.00011(0.788578889012330, -0.614933602760558)) C R 2 , 

and iterate each of them 10 times with F~ x . These 2200 points z^ n = F~ n (qi), 
1 < i < 200, < n < 10, provide an approximation of the stable manifold. 
We next construct rectangles 

Bi, n = {(x, u) E M 2 : \(x - V x T(zi, n ),u - V u T(z i>n )) • e"| < 0.001, 
\(x-V x T{zi, n ),u-V u T(z i>n )) -e s \ < 0.001}, 

where 

e u = (0.570868623900820281,-0.821041420541974730), 
e s = (0.992704972028756258,0.120568812341277774) 
are the approximate unstable and stable eigenvectors of DF* at pq, and verify that 

F F F F „ 

Si, io =^ Bi$ . . . ==>■ Bi o => So, 

for all F e F„, where 

B Q = {(x,u) Gl 2 : \(x~V x T( Po ),u-V u T(p Q ))-e u \ < 0.001, 
| (a; - r x T( Po ),u- V u T( Pa )) ■ e s \ < 0.001}. 

To prove that the cone conditions are satisfied we verify the hypothesis of Lemma 
using the routine checkConeConditions. This proves that ([J in Bi ^ U Bq 



6.11 



covers the local stable manifold. 

The sets ko,ko,kJ, and k*, from Lemma 11. l[ serve as bounds on the end points 
of the sets /Cg 1 in the Markov partition V = {/Cq , /C"}. Therefore, 

r L {V) > d ist(ky 4) ^ 0.0650166, (69) 

VI + Lip 2 sup gek t i<3 , ek b|g - q'\ 

r R {V) > . f st(k ° lM) > 0.0514139 (70) 

VI + Lip sup gek t i9 , ek b|£?- q'\ 

(these quantities are computed in the routine Gaps in (Progs2 2009)). 

Finally, we compute (in routine Distortion of (Progs2 2009)) the quantities C, 7 
and e appearing in the Definition 9.5 and verify parts 3), 4) and 5) of this definition. 
In particular, we have set 

f max{|6|,|c|H f , max{|6|, \c\} 

A UA! I \a\ - 1 J AqUAx I \a\ - 1 
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while 7 has been chosen as the supremum of all possible values that verify the 



conditions ( 52 1 — ( 53 1 over AqUA^ To verify part 4) we have bound the total variation 
of r](x, u) — log \a(x, u) | on A^ from above by Area(AQ • sup^^g A / | Vr/(x, u) \ . 

Lemma 11.2 Any function G £ G* : Ao U Ai i— > IR 2 is of class J~(C, 7, e) with 

C=16.6, 7 = 47.8, e = 0.88, 

A straightforward implementation of the Duarte Distortion Theorem and Lemma 
|9. 8 1 gives us a bound on the Hausdorff dimension of the unstable Cantor set: 

Corollary 11.3 For all G G G*, the distortion of the unstable Cantor set IC G is less 
or equal to 

L»(C,7,e) = 3749.3, 
while its Hausdorff dimension satisfies: 



l 0g (l I e- D rn(V) \ l 0g (l I e- D rL (V) \ 
iOg ^1 + 1+e o TL{v) J 10g ^1 + 1+e D TR(v) J 



1 „-2 D Tr{V) 



where tl(V) and tr are as in (69 1 -(70 1 



The Hausdorff dimension of the set Cq satisfies 

Aim H {C G ) > 2dim H (/C£) pb 0.00013 e" 7499 . 
Proof. The fact that the stable and unstable foliations are transversal (see Lemma 



10.41 implies that Cq = x ICq. Reversibility and the symmetry of the Markov 
partition Ao H Ai implies that JC G = T(IC G ), and dim//(/Cg) = dim#(/Cg.), which in 
turn implies the last claim. 

□ 

Corollary |1 1 .3 together with Theorem 10.3 proves the first part of the second 



Main Theorem. 

Let us now construct a convergent sequence of approximations of the Cantor sets 
C G . As before, we denote A and A^ the two full components of the Markov partition 
as in Definition |9.5| and Lemma [10. 5 [ A = Ao U Ai. Define recursively: 

U G = G(A)C)G~ 1 (A) and U G = Gipl^ 1 ) n G _1 (Wq _1 ). 



Each of the sets U G contains 2 • 4 components U G n , n = 1..2 • 4 
Recall that dAj = bj U tj U lj U Tj, j = 0, 1. By definition of U G ' n , the boundary 



dU G ' n consists of two images, t k ' n and b k ' n , of some subsets of tj and bj, j — 0, 1, 
under G k , and two images, l k ' n and r fc '™, of subsets of lj and rj, j = 0, 1, under G~ k . 
We will refer to t k ' n , b k,n , l k ' n and r k - n as "edges" or "sides" , and to their intersections 
as corners of U G . 

Lemma 11.4 Let 

and set pk = min ra {pfe in }. There exist constants C > and c > such that 



diam(z4' n ) <CAZ k = CK k 



Pk >cA k = 



+> 



where A± are as in (65) 
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Proof. The length of an edge of U G ,n is bounded by cA~ k for some A 6 A. This 
follows from the fact that Aj, j = 0,1, are constructed so that the tangent vectors to 
tj a nd bj, j = 0,1, are contained in C u , and those to lj and rj, j = 0, 1, — in C s (see 



(67 1 and (68 1), and from the fact that an edge of U G n is contained in the image of an 



edge of U G ' 3 for some j. 

Consider an edge olU G , say t k ' n . Set 

af = sup {maximum Zv : v G C u (p)} , = inf {minimum Zv : v G C u (p)} , 

where Zv signifies the angle between vector v and the x-axis, measured, say, 
counterclockwise. Similarly for b k ' n , l k ' n and r k ' n . Straightforward geometric 
considerations demonstrate that the angles af b l specify two quadrilaterals, 

containing and contained in Uq* 1 , respectively. 

Since the angle between any v G C u (p) and any u G C s (p) is bounded from below 



(see Lemma 10.41, the angles of the quadrilaterals are bounded away from and it, 
and their diameters satisfy: 

cA7_ k < diam (Inner Quadrilateral) , 

C A_ > diam (Outer Quadrilateral) . 

The claim follows. 

□ 

Corollary 11.5 The Hausdorff dimension of the Cantor set C G satisfies 

dim H (C G ) < bg( - 4 ) < 0.7673. 

- log K+ 

Proof. Clearly Cq C U g , for all k. As above, U G = Un=i^G™' wnere 

diam(Wg'") < const n\, for all 1 < n < 2 ■ 4 fc . 
Hence, for any k > 1, C%{C G ) < const ^nf, and if An d + < 1, then C%(C G ) = 0. 

□ 

12. Overview of the programs 

We will now give a brief overview of the programs used in the proofs. 

A subset of the domain of a map F G is determined using the program 



f indDomain of (Progsl 2009), used in Lemma 4.2 f indDomain discretizes the domain 
of the generating function in the x-domain, and computes rigorous bounds on the 
maximum and minimum of the generating function, s(y, x), on this slice of the domain. 
During the process of computing the maximum and minimum, the program also verifies 
that the conditions of the implicit function theorem are verified on the slice. 

To verify that a pair of h-sets satisfies the hypothesis of Theorem |6.6[ we use the 
routine checkCoveringRelations of (Progsl 2009). It discretizes the boundary of 
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a h-set, and verifies that each piece of the boundary is mapped so that the covering 
relations hold. 

To verify that a pair of h-sets with cones satisfies the hypothesis of Lemma f6.11| we 
use the routine checkConeConditions of (Progs 1 2009). It computes the derivative 
of the map on the initial h-set in appropriate coordinates, and checks the positive 
definite condition of the quadratic form V in the statement of the lemma using the 
Sylvester criterion. 

The sequences of h-sets used to enclose the stable manifold in Section [TT] are 
constructed in the program encloseStableManif old of (Progsl 2009). This program 
uses the routines checkCoveringRelations and checkConeConditions, to verify that 
on each sequence, the hypotheses of Theorem 6.6 and Lemma 6.11 respectively, are 
satisfied. 

As we have already mentioned parts of Theorem |10.3[ Lemmas |10.4| and |11.2| are 
proved with the help of the routine Distortion of (Progs2 2009). The input of the 
routine is a set of parameters that define the components Aq and A^ (see Table |3| and 
the cone fields (see ( 62 1— ( 63 1 ) . The routine verifies that the cone fields are invariant 
and computes the expansion rates ( 65 ) . It also computes the quantities C, 7 and e 
appearing in Lemma |ll.2| 

The thicknesses (|69[) and (70) of the Markov partition for the horseshoe of G* 
are computed in the routine Gaps of (Progs2 2009). 



13. Concluding remarks 

An obvious question is whether the Hausdorff dimension of the set Cq for all G in 
some subset of G* is independent of G. 

In the satellite work (Gaidashev and Johnson 2009) we were able to use a 
renormalization approach to show that the Hausdorff dimension of the "stable" set, 
that is the set on which the maximal Lyapunov exponent is zero, is indeed invariant for 
a subset of infinitely renormalizable maps. An essential ingredient of that proof is the 
fact that the renormalizations of all infinitely renormalizable maps in a neighbourhood 
of the fixed point converge to that fixed point on their domain of analyticity T>. 

In a somewhat similar fashion, the invariance of the Hausdorff dimension of the 
set Cq could be demonstrated if one can show that renormalizations of the infinitely 
renormalizable maps converge to G* on neighbourhoods Vk of the rescalings A^ g (Cg) 
of the hyperbolic sets. However, T> U / t = i.. 00 is not a connected domain, and 
convergence on T> does not imply that on Ufc = i..ooVfc. 
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